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Abstract. Let v be a real polynomial of even degree, and let p be the equilibrium probabil- 
ity measure for v with support S; so that, v(x) > 2 J log \x— y\ p(dy)+C v for some constant 
C v with equality on S. Then S is the union of finitely many bounded intervals with end- 
points 5j, and p is given by an algebraic weight w(x) on S. The system of orthogonal 
polynomials for w gives rise to the Magnus-Schlesinger differential equations. This paper 
identifies the r function of this system with the Hankel determinant det[J x^ +k p(dx)] 1 J 
of p. The solutions of the Magnus-Schlesinger equations are realised by a linear system, 
which is used to compute the tau function in terms of a Gelfand-Levitan equation. The 
tau function is associated with a potential q and a scattering problem for the Schrodinger 
operator with potential q. For some algebro-geometric q, the paper solves the scattering 
problem in terms of linear systems. The theory extends naturally to elliptic curves and 
resolves the case where S has exactly two intervals. 

MSC (2000) classification: 60B20 (37K15) 
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1. Introduction 

This paper concerns systems of orthogonal polynomials that arise in random matrix theory, 
specifically in the theory of the generalized unitary ensemble [26], and may be described 
in terms of electrostatics. We consider a unit of charge to be distributed along an infinite 
conducting wire in the presence of an electrical field. The field is represented by a real 
polynomial v(x) = ^2j =0 ajX J such that a2N > 0, while the charge is represented by a 
Radon probability measure on the real line. 

Boutet de Monvel et al [5, 29] prove the existence of the equilibrium distribution p 
that minimises the electrostatic energy. Under general conditions which include the above 
v, they prove that there exists a constant C v such that 

v{x) > 2 / log \x - y\ p{dy) + C v (x 6 R) (1.1) 
Js 
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and that equality holds if and only if x belongs to a compact set S. Furthermore, there 
exists g > and 

-oo < S 1 < 5 2 < S 3 < ... < S 2g+2 < oo (1.2) 



such that 



S = U 9 +l[S2j-i,S 2j ] (1.3) 



It is a tricky problem, to find S for a given and [10, Theorem 1.46 and p. 408] contains 
some significant results including the bound g + l<N + lon the number of intervals. 
When v is convex, a relatively simple argument shows that g = 0, so there is a single 
interval [21]. 

Definition. The n th order Hankel determinant for p is 



D„ = det 



/ x J+k p{dx) . (1.4) 



In section 3 we introduce the system of orthogonal polynomials for p, and in section 4, 
we regard _D n as a function of 5 = (5i, . . . , S2 g +2), and derive a system for differential 
equations for logD n , known as Schlesinger's equations. Let A(z) be a proper rational 
2x2 matrix function with simple poles at 5j; let aj be the residue at dj, and suppose 
that the eigenvalues of otj are distinct modulo the integers for j = 1, . . . , M. Consider the 
differential equation 

±* = A{z)*{z), (1.5) 

and introduce the 1-form 



1 M 

0,(8) = - ^ trace residue(A(z) 2 : z = Sj) d5j (1.6) 
to describe its deformations. Then O turns out to be closed by [19]. 

Definition, (i) The tau function of the deformation equations associated with (1.5) is 
t : C M \ {diagonals} — > C such that dlogr = O. 

(ii) Given a self-adjoint and trace-class operator K : L 2 (p) — > L 2 (p) such that < 
K < I, and P(t,<x) l ne orthogonal projection / i— >■ I[t j00 )/, the tau function of K is 
r(£) = det(I — P[ tj(X )K). (The definitions in (i) and (ii) are reconciled by Proposition 3.3.) 

In section 5 use the results from preliminary sections to prove that D n gives the 
appropriate r function for Schlesinger's equations. As an illustration which is of importance 
in random matrix theory, we calculate the tau function explicitly when p is the semicircular 
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law. When S is the union of two intervals, the Schlesinger equations reduce to the Painleve 
VI equation, as we discuss on section 7. Okamoto derived r functions for other Painleve 
equations in [27]. See also [15] and [2]. 

In sections 2,3 and 4, we develop standard arguments, then our analysis follows that 
of Chen and Its [8] , who considered the p that is analogous to the Chebyshev distribution 
on multiple intervals. Chen and Its found their tau function explicitly in terms of theta 
functions on a hyperelliptic Riemann surface; in this paper, we show by means of scattering 
theory why such solutions emerge. 

In terms of definition (ii), we have a real t(x) which is associated with a compactly 
supported real potential q(x) = — 2^-logr(2x) and hence the Schrodinger differential 
operator — + q(x). One associates with each smooth q a scattering function 0; then 
one analyses the spectral data of — + q(x) in terms of </>, with a view to recovering q. 
The Gelfand-Levitan integral equation links <p with q. 

In random matrix theory, tau functions are introduced alongside integrable kernels 
that describe the distribution of eigenvalues of random matrices, especially the generalized 
unitary ensemble; see [14, 26, 31]. Let X be a n x n complex Hermitian matrix, let A = 
(Ai < A2 < • • • < A n ) be the corresponding eigenvalues, listed according to multiplicity, and 
consider the potential V(X) = n~ x Y^=i v (^j)- Now let dX be the product of Lebesgue 
measure on the entries that are on or above the leading diagonal of X; then there exists 
< Z n < 00 such that 

^\dX) = Z' 1 e^{-n 2 V{X))dX (1.7) 

defines a probability measure on the nxn complex Hermitian matrices. There is a natural 

action of the unitary group U(n) on M n given by (U,X) h» UXU\ which leaves 
(2) 

invariant. Hence v\ is the generalized unitary ensemble with potential v. 

There exists a constant ( n such that p n (dx) = C~ 1 e~ nv ^ dx defines a probability 
measure on R; then we let Efc 1 be the orthogonal projection onto span{x J : j = 0, . . . , k— 1} 
in L 2 (p n ). The eigenvalue distribution satisfies 

J k{j : XjiX) < t}v<p(dX) = det(J - I (t ,oo)^") (1.8) 

where the right-hand side can be expressed in terms of Hankel determinants. For large 
n, most of the eigenvalues actually lie in S by results of [5, 26]. Moreover, there exists a 
trace-class operator K on L 2 (p) such that < K < I and 

det(J - I(t,oo)££")) -> det(J - l ( t,oo)*0 in ->• 00); (1.9) 
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call this limit r(t). Tracy and Widom [31] showed how to express such determinants in 
terms of systems of differential equations and Hankel operators. 
We introduce the matrix 

J 



-1 

1 



and apply a simple gauge transformation to (1.5). Then for a sequence of real symmetric 
2x2 matrices J/3 k (n), we consider solutions of the differential equation 

dx x-S k 

k=l 



and form the kernel 



Z(x) ^ (x ->■ Sj), 



K(x , y)= z Mim. (L11) 

y-x 



We show that the properties of K depend crucially upon the sequence of signatures of the 
matrices (Sj — Sk)J(3k(n)- In Theorem 8.3, we introduce a symbol function <fi from Z, a 
constant signature matrix a and a Hankel operator such that K = I^oT^. 

In section 9 we introduce <fi from (1.11), express <p m terms of a linear system as in [4] 
and hence obtain a matrix Hamiltonian H(x) such that 

r(2x) = exp(-f trace H(u) du\ (1.12) 

J X 



and prove that q(x) is meromorphic on a region. We regard —-r-z + q as integrable if 



i 2 

dip- 

—f" + qf = A/ can be solved by quadratures for typical A. This imposes severe restrictions 
upon q; indeed, Gelfand, Dikij and Its [12, 6] showed that the integrable cases arise from 
finite-dimensional Hamiltonian systems. In sections 10, 11 and 12 we consider cases in 
which — /" + qf = Xf has a meromorphic general solution for all A, and q satisfies one of 
the following conditions: 

(i) q is rational and bounded at infinity; 

(ii) q is of rational character on an elliptic curve; 

(iii) q is the restriction of an abelian function to a straight line in the Jacobian of a 
hyperelliptic Riemann surface. 

In cases (ii) and (iii), the corresponding Schrodinger equation has a spectrum with 
only finitely many gaps. In (i) and (ii), we introduce a linear system (—A,B,C) so as 
to realise <f>(x) = Ce~ xA B, and use the operators A, B, C to solve the Gelfand-Levitan 
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equation. Thus we obtain explicit expressions for and r. In (iii), we can do likewise 
under further hypotheses. 

2. The equilibrium measure 

Given the special form of the potential, the equilibrium measure and its support satisfy 
special properties. To describe these, we introduce the polynomial u of degree 2N — 2 by 

. . f v'(z) — v'(x) . , . 
u(z) = / —-^ —p(dx) (2.1) 

Js z ~ x 

and the Cauchy transform of p by 

f p(dx) {zeC \ s) (22) 
Js x - z 

and the weight 

23+2 



w(x) = 2Na 2N (-Q(x) JJ (x - S 2j -i)(x - <%)) ^ (2.3) 



where Q(x) is a product of monic irreducible quadratic factors such that w(x) 2 = 4u(x) — 
v'(x) 2 . 

Proposition 2.1 (i) The Cauchy transform is the algebraic function that satisfies 

R(z) 2 + v'(z)R(z) + u(z) = (2.4) 

and R(z) — > as z — > oo. There exist nonzero polynomials uo,iii,ii2 such that uqR' = 
uiR + u 2 . 

(ii) The support of p is 

S = {x E R: Au(x) -v'(x) 2 > 0}. (2.5) 

(iii) p is absolutely continuous and the Radon-Nikodym derivative satisfies 

where 2ir = f s w(t)dt and w(x) — > as x tends to an endpoint of S. 

Proof, (i) The quadratic equation is due to Bessis, Itzykson and Zuber, and is proved 
in the required form in [28]. One can easily deduce that R satisfies a first-order linear 
differential equation with polynomial coefficients. 
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(ii) Pastur [28] shows that the support is those real x such that 

\v'(x) + vV(» 2 - 4u(x)\ 2 = 4u(x), (2.7) 

and this condition reduces to 4u(x) > v'(x) 2 and u(x) > 0, where the former inequality 
implies the latter. The polynomial 4u(x) — v'(x) 2 has real zeros 5\, . . . , $23+2, and may 
additionally have pairs of complex conjugate roots, which we list as #2 g +3, . . . , 64^-2 with 
regard to multiplicity. Hence we can introduce w as above such that 4u(x)—v'(x) 2 = w(x) 2 . 
(hi) From (i) we deduce that 

since both sides are holomorphic on C \ S, vanish at infinity and have the same jump 
across S. By Plemelj's formula, we deduce that 



v (A) = 2p.v. - — — (XeS). (2.9) 



y/Mt) -V(t) 2 dt 

'S 

See [28, 5, 29]. This gives the required expression for p. 

□ 

3. Orthogonal polynomials 

First we introduce orthogonal polynomials for p, then the corresponding differential 
equations. Let (pj)^L be the sequence of monic orthogonal polynomials in L 2 (p), where 
Pj has degree j and let hj be the constants such that 

/ Pj(x)pk(x)p(dx) = hjd jk ; (3.1) 
Js 

and let (qj) ( jL 1 be the monic polynomials of the second kind, where 

w = f oW-aW ^, ( 3. 2 ) 

Js z-x 

has degree j — 1. On account of Proposition 2.1, the orthogonal polynomials are semi 
classical in Magnus's sense [22], although the weight typically lives on several intervals. 
The following result is standard in the theory of orthogonal polynomials; see [8]. 

Lemma 3.1 Let c n = h n /h n -i and b n = h~ x f s xp n (x) 2 p(dx). Then 

(i) the polynomials (p n )5^Lo sa tisfy the three-term recurrence relation 

Xp n (x) = p n +l(x) + b n+1 p n (x) + C n p n _i(x); (3.3) 
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(ii) the polynomials (qj)'jL 1 likewise satisfy (3.3); 
(Hi) the Hankel determinant of (1.4) satisfies 



Dn = ^0^1 • • • h n -\. 



(3.4) 



We introduce also 



Y n (z) = 



„ ( v \ C Pn(t)p(M) 

Pn\ z ) Js z-t 

Pn-l(z) 1 r p n -i(t)p(dt) 

h n —i h n —i JS 



z-t 



and 



V n (z) 



i/K o 



Proposition 3.2 (i) The matrices satisfy the recurrence relation 



(3.5) 



(3.6) 



Y n+1 (z) = V n (z)Y n (z). 



(3.7) 



(ii) The matrix Y n (z) is invertible, and dety n (z) = 1. 

Proof, (i) This follows from (i) and (ii) of the Lemma 3.1. 

(ii) This follows by induction, where the induction step follows from the recurrence 
relation in (i). 

□ 

We restrict p restrict to (— oo,t) H S and let 



(t) = [ x j p{dx) 



(3.8) 



be the j moment; the corresponding Hankel determinant is 



D n+1 (t) = det[p j+k (t)r k=Q . 



(3.9) 



Let E n : L 2 (p) — > span{x fc : k = 0, ...,n — 1} be the orthogonal projection; we also 
introduce the projection P(t, oo) on L 2 (p) given by multiplication / >->■ I(t )00 )/, where I(t j00 ) 
denotes the indicator function of (t, oo). 

Proposition 3.3 The tau function of E n+1 satisfies 



det(J- J B B+1 P (t , 00) ) = %^. 



(3.10) 



Proof. We introduce an upper triangular matrix [a>i,j]% =0 with ones on the leading 
diagonal such that Pj(x) = YH=o a tj xi - Then we can compute 



det[^+fc(*)]£ fc=0 = det[a £J ]l j=0 det[^ +fc (t)]" fc=0 det[a fc , m ]™ m=0 

n 

— oo ■*j>k=0 



= det 



Pj(x)p k (x)p(dx) 



(3.11) 



We can also express the operators on L 2 (p) as matrices with respect to the orthonormal 
basis (pj/ \/7i7)™ =0 , and we find 



E n +i — E n+ iP( tj00 -)E n+ i <H> 



1 



y/hjhh 



Pj(z)p k (z)p(dz) 



k J — oo 



j,k=0 



(3.12) 



so that 
det 

We deduce that 



/ Pj(x)p k (x)p(dx) 

J — oo 



j,k=0 



— det(E n+1 — E n+1 P( tyOO )E n+1 )h . . . h n . (3.13) 



det[p j+k (t)]™ k=0 = det(E n+1 - E n+1 P( ty00 ^E n+1 )D n+1 . 

4. Schlesinger's equations and recurrence relations 

Invoking Proposition 3.2(h), we introduce the matrix function 



(3.14) 



□ 



A n {z) =Y^(z)Y n (z)- 1 + Y n (z) 





-w'{z)/w{z) 



Yn(z)- 1 . 



(4.1) 



The basic properties of A n {z) are stated in (i) of the following Lemma, while (ii) gives 
detailed information that we need in the subsequent proof of Theorem 5.1. 

Lemma 4.1 (i) Let v'(z) 2 — Au(z) have zeros at 8j for j = 1, . . . , AN — 2. Then A n (z) is 
a proper rational function so that 



47V-2 



OLj{n) 



(4.2) 



where the residue matrices otj (n) depend implicitly upon S. 
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(ii) The (1, 2) and diagonal entries of the residue matrices satisfy 



4AT-2 



^2 a k( n )i2 = 0; 
k=l 

^2 ( a k( n )n - «fc(»22) = 2(n + N) - 1; 

4JV-2 

5fcafc(n)i2 = —2h n (n + N). 



4N-2 



k=l 



(4.3) 
(4.4) 
(4.5) 



k=i 



Proof, (i) The defining equation (4.1) for A n (z) may be written more explicitly as 



P'n{ z ) 

P'n-l( Z ) 



r p n (t)w(t)dt 
JS (z-t) 2 
l r p n -l(t)w(t)dt 
' h n -! JS (I^TP 



= A n (z) 



Pn(z) J, 
Pn-l(z) l r 

h n —i h n —i J i 



p n (t)w(t)dt 

S z-t 



p n - 1 (t)w(t)dt 
S z-t 



+ 



w'(z) r p n -i(t)w(t)dt 
w(z) JS z — t 



w'(z) f p n - 1 (t)w(t)dt 
h n -iw(z) JS z — t 



(4.6) 



By considering the entries, we see that A n (z) is a proper rational function with possible 
simple poles at the Sj, as in (4.2). Hence we have a Laurent expansion 



4N-2 



47V-2 



= - ak ^ + 4 5 kU k {n) + 0(4) (* ->• °o)- ( 4 -7) 
k=l k=l 



(ii) First we compute the (1, 2) entry of A n (z), namely 

4„(z)i 2 = -P' n (z) 

JS 



p n (t)w(t)dt _ f p n (t)w(t)dt _ /" p n (t)w(t)dt 

"z-t Pn[Z) .L {z-t)* ~w(z) Pn[Z) J s z-t 



z 



s * - L JS 
\ i t n p n (t)w(t)dt - in+ ^ 9 n(Z) I t"l>„(t)<r(t)<lt 



s 



j a 



w'(z) Pnjz) 

w(z) z n+1 



-p n (t)w{t)dt + o{^) 



(4.8) 



and we can reduce these terms to 



, , v -Kn h n (n+l) h n (2N-l) / 1 \ 



(4.9) 



which gives (4.3) and (4.5). 
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Next, the (2,2) entry of A n {z) is 

Pn-l( z ) f Pn(t)w(t)dt p n (z) f p n _i (t)w(t)dt 



A n (z) 



22 



h n -l Js z ~t h n-l Js ( z ~t) 

w'{z) p n (z) f p n - 1 {t)w{t)dt 
h n -i Js 

Pn(^J I _ ... | 



w(z) 
(n-l)z n - 2 



t n p n (t)w(t)dt - 



K_ lZ ^ J s w h n _ lZ n+iJ s 

2N-lp n (z) 1 



nt n ~ l p n -i(t)w(t)dt 



z z n h 

l-n-2N /l 
hOi 



- / t^Pn^W^dt+Of^ 
n-1 JS 



z \z* 
Similarly, the (1,1) entry is 



(4.10) 



PnW /" Pn-l{t)w(t)dt p n -i(z) f p n (t)w(t)dt 



z-t 



+ 



h n - X j s (z-ty 



+ 



W'{z)p n - 1 {z) f p n (t)w(t)dt 



w(z) h„ i J s : t 

Pn( Z ) 



[ t n - 1 p n _ 1 {t)w{t)dt + 
Js 



(n + l)p n _i(z) 



1 1„ 1 "-" Js ' h n l"-"' 2 Js 

W'(z) p n -l(z) 



t n p n (t)w(t)dt 



+ 



[ t n p n (t)w(t)dt + o(\ 
Js Kz 



w(z) /l n _i^ +1 Js 

= - z +0(^) (z->oo). (4.11) 
By comparing the coefficients of 1/z in (4.7) with (4.9), (4.10) and (4.11), we obtain 



4JV-2 



fc=i 



n 

1 - n - 2N 



which leads to (4.4). 



Let 



v / 27rip n (z) 

\/27Hp„-l(2) 



i-irw(z)p n (z)+q n (z) 
w{z)^/2iri 
i-Kw(z)p n - 1 (z)+q n - 1 (z) 
w(z)h n — i\/2iri 



(4.12) 



□ 



(4.13) 



which is a matrix function with entries in C(z)[n;]; note that <E> n also depends upon the Sj. 

Lemma 4.2 The functions $ n satisfy 
(i) the basic differential equation 

d$ n (z) 



dz 



= A n (z)$ n (z), 



(4.14) 
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(ii) the deformation equation 



<9$ r 
~86~ 



ctj (n) 



> 'J Z ~ 5 J 

(Hi) and the recurrence relation $ n+ \(z) = V n (z)3> n (z); 
(iv) moreover, $ n is invertible since det$ n (^) = l/w(z). 

Proof, (i) We can write 



(4.15) 



$ n (z) = Y n (z) 



V^d 
l/(w(z)V2wi) 



(4.16) 



and then the property (i) follows from (4.1). 

(ii) This follows from (i) by standard results in the theory of Fuchsian differential 
equations as in [14, 16]. 

(iii) The recurrence relation from Proposition 3.2(i). 

(iv) Given (iii), this identity follows from Proposition 3.2(h). 

□ 

Lemma 4.2 states several properties that the $ n satisfy simultaneously, and hence generates 
several consistency conditions. By taking (i), (ii) and (iii) pairwise, we obtain three Lax 
pairs, which we state in the following three propositions. 

Proposition 4.3 The residue matrices satisfy Schlesinger's equations 

da k (n) _ [aj(n),a k (n)} 



and 



86, 



dctj (n) 
~86~ 



So - 6k 



U * k) 



4JV-2 



- £ 



Sj - 6 k 



(4.17) 



(4.18) 



Proof. We can express the consistency condition 9 d ^"g^ = 9 as the Lax pair 



&A n (z) _ o^ra) _ OLj{n) aj{n)A n (z) 



86, 



Z — 6n 



(z - 5jY 



z-Sj 



(4.19) 



and then one can simplify the resulting system of differential equations. See [14, 19]. 



□ 



Proposition 4.4 The basic differential equation (4.14) and the recurrence relation in 
Lemma 4.2 are consistent, so 



A n+1 (z)V n (z) -V n (z)A n (z) 



1 




(4.20) 
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Proof. The Lax pair associated with the these conditions gives 

A n+1 (z)$ n+1 (z) = -^-$ n+1 {z) = -^-(y n {z)$ n (z)y 



(4.21) 



□ 



Proposition 4.5 (i) The deformation equation (4.15) and the recurrence relation in 
Lemma 4.2 are consistent, so 



aj{n+1) -V n (z) + V n (z) a > {n) - dVM 



z-Sj 

(ii) In particular, the (1,2) entry satisfies 

d 



Z — 8n 



85, 



85, 



logh n = -^ n 1 a J -(n)i 2 . 



Proof, (i) This is the Lax pair associated with Lemma 4.2. 
(ii) By letting z — > oo in (4.22), we deduce 

I~1 nl Ti nl T - dbn+1 

- aj (n + l) n n + n n «» = 

.9. as yj 



dSi 



K d& 3 



which implies that ctj (71)12 = 



dh r , 



5. The tau function 

We introduce the differential 1-form on Q 4N ~ 2 \ {diagonals} by 

4N-2 



trace 



V 5j -5k ' 



Theorem 5.1 The Hankel determinant D n gives the tau function, so 

Q n = dlogD n . 



(4.22) 



(4.23) 
□ 



(4.24) 



□ 



(5.1) 



(5.2) 



Proof. By Proposition 4.3 and results of Jimbo et al [20], fi n is an exact differential form, 
so dQ n = 0; hence there exists a function r n such that <ilogr n = Q n , and so we proceed 
to identify r n . By Lemma 3.1(iii), we have \ogh n = logD n+ i/D n , hence we consider 



4AT-2 



n n+1 -n n = V tracef ^ (n + 1) " fc(? ! + 1 ]-^ (n)afc(n) )^ (5.3) 

j^k:j,k=l 
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where by Proposition 4.5(i) ctj(n + 1) = V n (5j)aj(n)V n (5j) 1 so 

trace(aj(n + l)ak(n + 1) — aj(n)ak(n)) 

= trace(a i (n)y n (5 j )" 1 y n (4)a fc (n)l/ n (4)- 1 y n (5 j ) - ctj (n)a k (n)). (5.4) 

We have 

V n (5 3 )- l V n (5 k ) = \ Sj -s k J] (5.5) 

so by direct calculation 

4JV-2 

O n+ l-O n = |/ln 1 «j( n )l2(«fc(^)ll - «fc(^)22) 

+ h~ 1 a k (n) 12 (a j (n) 2 2 - aj(n)u) 

- h~ 2 (5j - <yfc)aj(n)i 2 afc(n)i 2 }d5j (5.6) 

In this sum we have taken j 7^ k, but the expression is unchanged if we include the 
corresponding terms for j = k; hence the coefficient of dSj is 

4N-2 4N-2 



«j(™)i2 ^2 }l n 1 ( a k( n )n - (71)22) - {aj(n)n - a •/ (71)22) ^ h nak{n)i2 
fc=i fe=i 

r / \ 47V-2 / s 4JV-2 

-^^2 2^ «fc( n )l2 + ^T2 2^ ^«fc( n )l2- ( 5 - 7 ) 

n fc=l n fc=l 

We use Lemma 4.1 to reduce this to —h^ 1 etj(n)\2, so 



4AT-2 



O n+ i - fi n = - /i n 1 a J -(n)i 2 d<Jj 



47V-2 



^ Al g/, n ^.. (5.8) 



i=i J 



Hence Vt n = d Y2=l lo § 



□ 



Following [19], we interpret (5.1) in terms of integrable systems and Hamiltonian 
mechanics. Let M = M 2 (R) 47V ~ 2 be the product space of matrices, and let G = GL 2 (H) 
act on M by conjugating each matrix in the list 

(X 1 ,...,X n )^ (UX&- 1 , UXnU' 1 ). 



13 



The Lie algebra g of G has dual g* , and for each (eg* the symplectic structure at £ on 
g x g is given by u^(X, Y) = £([X, Y]). Given 



2N-2 



A M = £ ^ ,5.9) 

fc=l 

as in (4.2), we introduce 

2N-2 

u(X,Y)= trace(a k [X k ,Y k }) (5.10) 



k=i 



for X = (X k ) 2 k ^- 2 and F = (Ifc)^ 2 in S 2N ~ 2 - Given f,g : M ^ C, their Poisson 
bracket is {/, = Xf{g), and the corresponding vector field satisfies Xr/, s } = LY/, 
The spectral curve of A(z) is the algebraic variety 

E A = w)eC 2 : det(wJ - = o}. (5.11) 

As suggested by (5.1), we introduce the Hamiltonian H 3 ■ : M — > C by 

Hj = E (, 12) 

so that ^7 logr(<5) = Hj. We observe that Hj is a polynomial in the entries of aj(n) and 
a k (n), and is a rational function of 5j and 5fc. To lighten the notation, we temporarily 
suppress the variable n. 

Proposition 5.2 (i) The Hamiltonian Hj gives a vector Geld (X^) 2 ,^ 2 which is associ- 
ated with the differential equation 



dA 
~dt 



z-Sjl 



(5.13) 



(ii) The Poisson brackets of the flows commute, so that {Hj, H k } = 0. 
(in) Under this How, the spectral curve of A is invariant. 

Proof, (i) For each Y = (Y k )™- 2 , we introduce a flow on M by a k = [Y k , a k }. We can 



differentiate Hj in the direction of Y and obtain 



Y(H )= V trace ([y fc ,a fc ]o!j) | trace (a k [Y k , Oj\) 
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With H k , we associate the Hamiltonian vector field X H . = (X { H ') k=1 such that 

27V-2 

Y(H j ) = u(X H .,Y) = tmce(a k [(X^,Y k }). (5.15) 



k=i 

We deduce that 

< = 5^ (6 ' 16) 

x «) = £ jp^- (5-17) 

It is then a simple calculation to check that a k = [X^,a k ] extends to give (5.13) for 
(d/dt)A(z). 

(ii) Given the vector fields (Xjp)j corresponding to H k and (Xjp)j corresponding to 
He from (5.12), one can compute 

{H kl H t } = ^trace([X^,^]xg) (5.18) 



and reduce the expression to zero by an elementary calculation. 

(iii) One can check that for each positive integer m, the ^traceA(z) m = 0, and hence 
det(wl — A(z)) is invariant under the flow. 

□ 

6. Orthogonal polynomials on a single interval 

In this section we consider the Chebyshev polynomials. 

• First suppose that v = 0. Then the corresponding equilibrium distribution is the 
Chebyshev distribution (l/-7r)(l — a; 2 ) -1 / 2 . In this case, the orthogonal polynomials are the 
Chebyshev polynomials of the first kind and, unsurprisingly, our results reduce to those of 
Chen and Its [8]. 

• For a < b, let 



(b-a) 2 \ 2 

so, by standard results used in random matrix theory [26], the equilibrium measure is the 
semicircular law on [a, b], as given by 

8 



p(dx) = — ^y/(b-x)(x-a) I[ ayb ] (x)dx (6.2) 

n(b — ay 
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Proposition 6.1 The tau function for the semicircular distribution is 

T(a, b) = (a- b )(^+2n+l)/4 e n(n + 2)(a-bf /32 (g^ 

Proof. Let t/ n be the Chebyshev polynomial of the second kind of degree n, which satisfies 

, _ . sin(n + 1)6> , 

Z7 n cosfl = . Q , 6.4 
sine* 

and let 

p n (x) = 2~ 2 "(b - a) n U n ( 2X ~i a a +5) ) (6-5) 

which is monic and of degree n, and the p n are orthogonal with respect to the measure 
p. By elementary calculations involving trigonometric functions, one can show that h n = 
2- 4n (b-a) 2n and 



A n (x) ' 



(x — b) (x — a) 



n—li 



n(x- (a + b)/2) -(n + 1)(6 - a) 2 h 
n{b - a) 2 /(2/i n _x) -\n + l)(x - (a + 6)/2) 



(6.6) 



which has poles at a and 6, as expected. One verifies that 



/n 2 + (n + 1) 2 \ da - db n(n + 2). „, . . 

^n=( ^ -) a _ b + 16 J (a - b)(da - db) (6.7) 

so that (6.3) follows by integration. 

7. Painleve equations for pairs of intervals 

Akhiezer considered a generalization of the Chebyshev polynomials to the pair of intervals 
[—1, a]U[/3, 1], and investigated their properties by conformal mapping. Chen and Lawrence 
[9] used the theory of elliptic functions to investigate these polynomials and in (8.18) 
expressed the Hankel determinant in terms of Jacobi's elliptic theta functions. In this 
section we obtain the differential equation where S is two intervals, and obtain a differential 
equation for the endpoints that is related to the one from [9] . 

Let v be a polynomial of degree 2N > 4 such that S = [Si, 82] U [£3, 84]. There exists 
a Mobius transformation ip such that (p(8i) = 0, ^(82) = 1 and ^(84) = 00; then we let 
t = <p(8s). Having fixed three of the endpoints, we can introduce the differential equations 
from section 4 that describe the effect of varying the endpoint t, namely 

d ^ . . fan ai «t \, /„„x 

and 

— = 7.2) 

dt x-t K 1 
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Let A(x,t) be the matrix (oto/x + ot\/{x — 1) + a t /(x — t)) and let A(x,t)i 2 be its top 
right entry. Then we introduce x = X(t) such that A(x,t)i 2 = 0; then by [18, p. 1333], 
the corresponding Schlesinger equations give a version of the nonlinear Painleve equation 
Pvi in terms of A, namely 

d 2 X (\ 1 1 \dX 1/1 1 1 \/dA\ 2 

+ -: r + 



eft 2 Vt £-1 X-tJ dt 2 VA A-l X-tJ\dtJ 

_ lA(A-l)(A-t) / M fei(*-l) (fct -!)*(*- 1)' 



a 00 a 2 + a-D 2 fA-ti 2 



2 t 2 (t-l) 2 V 00 A 2 (A-l) 2 (x-ty 

The Hamiltonian and tau function satisfy 

H t = - log r = trace (— + —). (7.4) 

Having transformed S 1 to [0, 1] U [t, 00] we can lift this to the portions of the real axis that 
are covered by the elliptic curve £ = {(A, w) : w 2 = 4A(A — 1)(A — t)} which has parameters 
A = V{u/2) and w = V'{u/2) in terms of Weierstrass's function V with e\ = t, e 2 = 1 and 
es = 0. Hence we transform to the dependent variables 



u= I - . "° w = , (7.5) 
Jo 



(is 

v/s(s^l)(s-t)' 



and make the substitution u = u(X(t),t). Fuchs [16] observed that 
d 2 u 2t — 1 du u 



dt 2 t(t - 1) dt 4t(t - 1) 

y/X(X-l)(X-t) \ k t k x {t - 1) k t t(t-l) 

•Zoo ~Tn r 



2t 2 (t-l) 2 L A 2 (A-l) 2 (A-t) : 
To solve this in a special case, we introduce the complete elliptic integral 



(7.6) 



if(m) = / 



a/1 — m 2 sin 2 

By comparing terms on power series, one can recover the following result. 

Proposition 7.1 (Poincare) Suppose that ko = k\ = k t = k^ = and that u(t) = 
ciK(\/t) + c 2 K{\/\ — t) for constants c\ and c 2 . Then u satisfies Legendre's equation 

, d 2 u . du u , . 

<(< - 1) ^ + (2i - 1 )* + 4=°- (7 ' 7 » 

so X(t) = V(u(t)/2; 0, 1, t) gives a solution of Pyi- 
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8. Kernels associated with Schlesinger's equations 

In this section, we introduce kernels that are associated with Schlesinger's equations, and 
then factorize them in terms of Hankel operators. First we let Vj = — 2 _1 trace ctj(n) 
and observe that Vj does not depend upon n. Indeed, by multiplying (4.22) by V~ x , 
one deduces that trace A n+ i(z) = trace A n (z), and since trace ctj(n) = \im z ^.s ( z — 
5j) trace A n {z), we deduce that trace ctj(n) is constant with respect to n. By (4.12), we 
have Xlj=i~ 2 trace aj(n) = 1 — 2N. Now, given $ n as in (4.13), let 

4N-2 

*»(*) = J] (Z-S^^Z). (8.1) 
3 = 1 

We next introduce the matrix valued kernel 

M ^' C) - -27ri(z-C) ' ( } 

we aim to show that M n is positive definite as an integral operator on L > (S), and we 
observe that this property does not change if we introduce weights on S. 

Proposition 8.1 Let E n (z,Q be the kernel of the orthogonal projection onto span{x J : 
j = 0, . . . , n — 1} in L 2 (p). Then the top left entry of M n (z, Q equals 

, 4N-2 4N-2 

M„(z,C)n = t-=- I] (*-*i) Vi II (C-SjpE n (z,0. (8.3) 

Aln " 1 j=l 3=1 

Proof. The Christoffel-Darboux formula gives 

tp fv *S Pn(z)Pn-l(Q - Pn-l(z) Pn (() 

En{z > = hJz—Q ' ( } 



One can find ^ n (z)^ J^ n (() by direct calculation, and compare with this. 



□ 



Let (3j(n) = osj(n) + Vjlii which has zero trace. Furthermore, if $ n is a solution of 
the basic differential equation (4.14), then 

jfn{z) = B n (z)V n (z) (8.5) 



where 

4N-2 



B »M = E 7=X- (8-6) 



Z — 5n 

3 = 1 3 
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We pause to note an existence result for solutions of the matrix system (8.5). 

Lemma 8.2 Suppose that (3j(n) has eigenvalues ±Kj(n) where 2nj(n) is not an integer. 
Then on a neighbourhood of Sj, there exists an analytic matrix function S n j such that 



* n (z) = E n , j (z)(z-6 j )l 

satisfies (8.5). 

Proof. This follows from Turrittin's theorem; see [3]. 



(n) 



(8.7) 



□ 



For notational simplicity, we consider the interval (Si, 62) and assume that 5\ = and 
1 < 5 2 ; the general case follows by scaling. For a continuous function : (0, 1) — > R 87V ~ 6 , 
the Hankel operator r> : L 2 ((0, l);dy/y; R) ->■ L 2 ((0, 1); dy/y; K 8N ~ 6 ) is given by 



T^f(x) = [ <j>(xy)f(y) 
Jo 



y ' 



(8.8) 



Since (3k( n ) has zero trace, the matrix (6\ — 5k)J/3k(n) is real symmetric and hence is 
congruent to either 



1 
1 



± 



1 




1 
-1 








(8.9) 



let a = diagonal [u/s] ^ 2 2 be the block diagonal sum of these matrices. 

Theorem 8.3 (i) Let /?i(n) be as in Lemma 8.2. Then there exists Z n , a 2 x 1 real vector 
solution of (8.5) such that Z n (x) — > as x — > 81 . 

(ii) The integral operator on L 2 ((0, 1); dx/x) with kernel 



K n (zX) 



(Z n (()UZ n (z) 
z-C 



(8.10) 



is of trace class; moreover, there exists a real vector Hankel operator T^ n on 
L 2 ((0,l),dy/y) such that 



K„ 



(8.11) 



(in) lfa>0, then K n >0. 
Proof, (i) There exists an invertible constant 2x2 matrix S n such that 



SnZ^S- 1 



z Ki(n) q 
Z -Ki(n) 



(8.12) 
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where Ki(n) > 0. Hence by Lemma 8.2, there exists a constant 2x1 matrix C such that 
Z n (z) = ^ n (z)C is a solution of (8.5), and Z n (z) = 0(\z - d^ 1 ^) as z ->■ 5i. 

(ii) Hence we can introduce K n by (8.10), and next we prove that the kernel satisfies 



First note that by homogeneity (x-^ + y-^)^/xy/(x — y) = 0. Since the fikip) have zero 
trace, we have Jf3k(n) + (3k(n)^ J = and hence the differential equation gives 

(x— + y—^jZ n {y)^JZ n {x) 



= Z n (y^B n {yfjZ n (x) + Z n (y)^ JB n {x)Z n {x) 

' x y 

^x - 5 k y-5 k 



4N-2 



^yz n ( y yj^(n)z n (x)( v ^ J - - -^-); 



k=2 



(8.14) 



on dividing by x — y, we obtain 

d d \y/xyZ n (y)^JZ n (x) ' A 



(6 6 \ 
\ X d-x +V dy) 



x-y 



= E 



hJxy 



k=2 



(x - S k )(y - S k ) 



Z n (y)^jp k (n)Z n (x) (8.15) 



as in (8.13). 

Noting the shape of the final factor in (8.12), we choose 



'Pni.x) = column 



y/xZ n {x) 

x — 8h 



k=2,...AN-2 



(8.16) 



which has a 2 x 1 entry for each endpoint 5 k of S after 8i, and the block diagonal matrix 



P(n) = diagonal —S k J Pk(n) 
with 2x2 blocks, and we consider 



fc=2,...,4JV-2 



dz 



(8.17) 



(8.18) 



K n (x,y)= (j) n (yz)^p(n)(j) n (zx) — . 

Jo z 

First note that since K\{n) > 0, we have K(x, y) — > as x, y — > 0. Then 

i X dx +y dy)^ n(yX ' y ^ = J ( y 0'n( yz ^P( n ^ n ( zx } + X( Pn(yzy (3(n)(j)' n (xy)) dz 

= <j> n {y)^{n)<j> n {x) - 0„(O)+/3(n)0 n (O). (8.19) 
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We have n (O) = 0, so 

K n {x, y) = K n {x, y) + £(x/y) (8.20) 

for some function £. But Z n (z)/(z — 8\) Kl ^ n ' is analytic on a neighbourhood of <5i, so it is 
clear that K n (x, y) — > and K n (x, y) — > as x — >■ or y — > 0; hence £ = 0. 

By the choice of ex, there exists a block diagonal matrix 7(n) such that 7(n) 'o"7(n) = 
/3(n), so we can introduce ip n (x) = 7(n)</> n (x) such that (f) n (xy (3(n)(f) n (y) = ip n (xY aip n (y). 
For this symbol function tfj n we have 

f 1 dz 
K n (x,y)= ifj n (yz) j atfj n (zx) — , (8.21) 
Jo z 

or in terms of Hankel operators K n = crT^ n . We have 

/ \og\u n (u)\\ 2< ^- < 00, 

so r.^ is Hilbert-Schmidt and hence K n is trace class. 

(iii) If crfc > for all or equivalently a > 0, then K n > 0. 



□ 



Corollary 8.4 Suppose that Z is a 2 x 1 solution of 



±Z=(^ + ^ + ^-)z (8.22) 
dx \ x x — 1 x — t) 

such that the entries satisfy Z(x) = Z{x) and where 

(i) Po is as in Lemma 8.2, and Z(x) — > as x — > 0; 

(ii) J /3i is positive definite; 

(iii) jp t > 0. 

Then there exist an invertible real matrix S and a real diagonal matrix D such that 



^/xSZ(x) 

y/iD5Z(x) 
x—t 



(8.23) 



satisfies ip(x) = ij)(x) and 

JxyZixjf J Z{x) 



y-x 



f 1 dy 

= / ip(xz)^(zy) — . (8.24) 
Jo z 



Proof. We simultaneously reduce the quadratic forms associated with J Pi and J(3 t , and 
introduce an invertible real matrix S such that J Pi = SS^ and J fit = SD 2 S^ where D 
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is a real diagonal matrix such that the diagonal entries k of D 2 satisfy det(/3i — K/3 t ) = 0. 
Then we can write 

J Pi 



^(y)V(x) = y/xyZiv)* ( j^—J^ZTY) + {x - t){y - t) ) Z {x) ' 
Now we can follow the proof of Theorem 8.3, and deduce that 



(8.25) 



-i/>(y)ty(x) = fx— + y 



d d \ y/xyZ(y)UZ(x) 

' k 7/ — • 



dx dyJ x — y 
hence we can obtain the result by integrating and using (i). 

Example 8.5 For the semicircle law on [a, 6], as in Proposition 6.1, we have 

i(b-a) 2n+l 



and 



so that 



Jf3 a (n) 



Jf3b(n) = 



2/i n _i 4 
2n+l (n+l)(b-q)fe ra _i 
4 8 



-n(b— a) 

2/ln-l 

2n+l 



2n+l 
4 

-(n+l)(b— q)/t ra _i 



det J(3 a (n) = det J/?b(n) = 



n(n + !)(&- a) 2 (2n + 1) ; 



?.26) 



□ 



(8.27) 



(8.28) 



1.29) 
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In particular, when a = — 1 and 6=1, the matrices Jj3-i(n) and Jj3i(n) are indefinite. 
9. The tau function realised by a linear system 

In this section, we express the tau function of K n from Theorem 8.3 as a Fredholm 
determinant, and then obtain this from the solution of an integral equation of Gelfand- 
Levitan type. The first step is introduce a scattering function ip and then to realise this 
by a linear system, so that we can solve the Gelfand-Levitan equation. 

The differential equation 



dx 



B n {x)Z n {x) 



has a solution from which we constructed a symbol function 

"v^7(n)Z(x)n4iv-2 



V'n(^) = column 



k=2 



(9.1) 



(9.2) 



x - 5 k 

Suppressing n for simplicity, we change x G (0, 1) to t e (0, oo) by letting x = 8\ +e~ f and 
in the new variables write 

oo 

m = J2xie~^ +e+1/2)t . (9.3) 
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where YltLo \\Xi\\ < 00 • Likewise, we write r(t) for t(Si + e *). 
Let = {z : Ute > 0} be the open right half-plane, let 



^(x) 




(9.4) 



and extend \& to an analytic function \& : fi — >■ Mgjv-5(C) such that \I/(:r) = ^(xY for 
x > 0. Let ^( s ) = ty(x + 2s) and ^*^(x) = ty(x + 2s) t and let a be a constant matrix; 
then let K s = F^* oT^ ( . be a family of operators on L 2 (0, 00). 

Proposition 9.1 (i) The r function associated with K = r^»oT^ is r(2s) = det(I — K 8 ), 

which gives an analytic function on O. 

,2 

(ii) Let q(s) = —2-^ logr(2s). Then q(s) is meromorphic on O, and analytic where 
x\\y(x + s)\\ 2 dx < 1. 

(in) If < K < I, then r(s) is non-negative for < s < 00, increasing and converges 
to one as s — > 00. 

Proof, (i) The kernel of the Hankel operator Tq,, has a nuclear expansion 



F *( ) ^2e~ {Kl+i+1/2){x+y+2s) 



£=0 



Xi 

x\ 



(9.5) 



w 



here YlT=o llX-dl J Q °° e _2( -' tl+ ^ +1 / 2 - )< ' x+3ffs - ) (ia; < 00, so the Fredholm determinants are well 
defined. As in Schwarz's reflection principle, s 1— > Vl/j^ is analytic, and r^, (s) is Hilbert- 
Schmidt, so K s is an analytic trace-class valued function on O. Using unitary equivalence, 
one checks that 

det(J - K a ) = det(J - P (2s ^K) (s > 0). (9.6) 

(ii) Except on the discrete set of zeros of r(2s), the operator / — K s is invertible and 

rf.) = ,9.7) 

(iii) This follows from (9.6). 

□ 

Next, we obtain an alternative formula for q by realising \1/ via a linear system. The 
technique is suggested by the inverse scattering transform. Let H = C 8N ~ 6 be the column 
vectors, H = £ 2 be Hilbert sequence space, written as infinite columns, and introduce an 
infinite row of column vectors C G £ 2 (H ) by C = (xt/WxeV^^iLo an d a column Bel 2 
by B = (\\xt\\ 1 ^ 2 )eLo an d the infinite square matrix A = diagonal [£+ K\ + 1/2]£L . While 
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A is real and symmetric, we will write in some subsequent formulas, so as to emphasize 
their symmetry. 

In the following result we use the (8N — 5) x (8N — 5) block matrices 



U(x,y) v(x,y) 
w(x,y) ] z(x,y) 



W(x,y) = 

so that ^(x) = ^(x)^ and the matrix Hamiltonian 

H(x) = 



ip(x) 
ifj(x)^ 



(9.8) 



U(x,x)a v(x,x) 
w(x,x)^(j z(x,x) 



(9.9) 



" U v~ 




' 






r vip* 




z 


* 














where v, w G Hq, U operates upon Hq and z is a scalar. To simplify the statements of 
results, we use a special non-associative product *, involving a, that is defined by 

(9.10) 

Theorem 9.2 (i) The symbol ip is realised by the linear system (—A, B, C), so 

</,(£) = Ce~ tA B. (9.11) 
(ii) There exists a solution of the Gelfand-Levitan equation 

poo 

W(x,y) + ^(x + y) + W(x,s)*^(s + y) ds = (0 < x < y) (9.12) 

J X 

such that the tau function of Proposition 9. 1 (i) satisfies 

d 



dx 



logr(2a;) = trace H(x) (x > 0). 



(9.13) 



(in) Suppose moreover that x\\V(x)\\ 2 dx < 1. Then 



(9.14) 



Proof, (i) This identity follows from (9.3). Since + £ + 1/2 > 0, the semigroup 
e~ tA = diagonal [e-H^i+^+V 2 )]^ Q consists of trace class operators, and the integrals in 
the remainder of the proof are convergent. 

(ii) We introduce the observability Gramian 



Q 



POO 

1=1 e~ sAt C^aCe~ sA ds (x > 0), 

J X 



(9.15) 
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modified to take account of ex, and the usual controllability Gramian 



J X 



e - sA BB ] e 



ds, 



(9.16) 



both of which define trace class operators on £ 2 , and where L x > 0. The controllability 
operator S x : £ 2 ((0, oo); Hq) — > H is 



/•oo 

E x f = / e- tA Bf(s)ds 

J X 

while the observability operator is 6 X : L 2 ((0, oo); H ) — > H is 

/•OO 

e K /= / e - sAt cV(s)rf s . 



(9.17) 



(9.18) 



Finally, we let ip^(s) = ip(s + 2x), so that ^(x) is realised by (— A,e xA B,Ce xA ). In 
terms of these operators, we have the basic identities 



t _^t e 



while 



-^x — "X"x 



and Q£ = Q x a&l. 
Hence we can rearrange the factors in the Fredholm determinants 



(9.19) 
(9.20) 



We deduce that 



det(J - XT^aT^) = det(J - ASt Q x aQ x E x ) 

= det(J - XE x ElQ x aQl) 
= det(I-XL x Ql). 



\ogr(2x) = logdet(J- r^ar^P (2Xi0o) ) 
= log det(J - ar^P (2Xi0o) r^) 
= logdet(/-ar^ ) r^ ) ) 

= io g det(/-r; M ar^ ) ) 

= trace log(J - L x Q a x ), 



(9.21) 



(9.22) 



and hence 
d 



^ logr(2x) = trace ((/ - L x Q a x )- x {e~ xA BB^ e~ xA Q a x + L x e~ xA " C ] aCe~ xA )^ 



fit, 



-xA 1 " 



QT x {I-L x Qr x )- x e-* A B 



+ trace aCe~ xA (I - L x Ql)- x L x e- xA * CK 

25 



(9.23) 



The integral equation 



+ 



+ 



U(x,y) v(x,y) 
w(x,y) ] z(x,y) 

U(x,s) v(x,s) 
w(x,s)^ z(x,s) 



ip(x + y) 

ip(x + y)^ 

i>(s + y) 
ifi(s + yy 



ds = 



reduces to the identities 

poo 

U(x,y) = - v(x,s)ijj(s + y) Jf ds, 

J X 

poo 

z(x,y) = — / w(x, s) ^aip(s + y) ds, 

J X 

and the pair of integral equations 

poo poo 

v(x, y) + ip(x + y) — / / v(x,t)ip(t + s)^ aip(s + y) dsdt = 

J X J X 

and 



oo poo 



w(x, y) + ip(x + y) — j \ ip(s + y)ip(t + s)^aw(x,t) dtds = 0. 

J X J X 

To solve these integral equations, we let 

v(x, y) = -Ce~ xA (I - L x Q° x y\-y A B 



and 



w 



(x,y) = -Ce-y A {I-L x QlY x e- xA B- 



(9.24) 



(9.25) 



(9.26) 



(9.27) 



(9.28) 



(9.29) 



then by substituting these into (9.25) we obtain the diagonal blocks of the solution W, 
namely 

U(x, y) = Ce~ xA (I - L x Q° x y x L x e-y A " (9.30) 



and 

z(x, y) = B^e-y A 'Ql{I - ^QD^e'^B. 
Hence we can identify the trace of the solution (9.9) as 

trace H(x) = trace all (x, x) + z(x, x) 

= tvaceaCe- xA (I - L x Q x )~ 1 e~ xAt 
+ B^ e~ xAi Q X (I - L x Ql)- l e~ xA B 

= ^logr(2x). 



(9.31) 



(9.32) 
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(iii) By integrating by parts, we obtain the identity 

for < x < y. One can easily verify that the product * and the standard matrix multipli- 
cation satisfy (QW) * \& = QiW * \&), hence the formula 

-2^-W(x,y) -2^y(x + y) - [ (2^-Wlx, sj) * *(s + y) ds = (9.34) 
dx dx J x V dx / 

follows from multiplying (9.33) by —2^, and this shows that both —2^W{x,y) and 
— ^')^( a; ' 2/) are solutions of the same integral equation. By uniqueness of solutions, 
they are equal. 

□ 

10. Integrability of the tau function of a linear system 

Let q(x) = — 2-i^ traced (x) and r be as in (9.32). In this section we describe the 
properties of r in terms of the algebraic theory of differential equations [30]. 

Let F be a field (of complex functions) with differential d that contains the subfield 
C of constants and adjoin an element h to form F(h), where either: 

(i) h = f g for some g G F, so dh = g; 

(ii) h = exp J g for some g G F; or 

(iii) h is algebraic over F. 

Definition. Let Fj (j = 1, . . . , n) be fields with differential d that contain the subfield C 
of constants and suppose that 

Fi CF 2 C...CF n , (10.1) 

where Fj arises from Fj_i by applying some operation (i), (ii) or (iii). Then F n is a 
Liouvillian extension of Fi. 

Example. The tau function (6.3) belongs to some Liouvillian extension of C(x). 

Lemma 10.1 Let A n : C n ->■ C n , B n : C m ->■ C n and C n : C n ->■ C k be finite matrices, 
such that MA n + A\M is positive definite for some positive definite M. Let ip n (t) = 
Ce~ tAn B, and define the corresponding terms as in the proof of Theorem 9.2. Then r n (x) 
belongs to some Liouvillian extension field F of Fo = C(t, e~ tK: > , e~ tKj ;j = 1, . . . , n), where 
( K j)j=i is the spectrum of A n . 

Proof. By Lyapunov's criterion [7], all of the eigenvalues Kj of A n satisfy JR/tj > 0, hence 
|| e -tA„|| j g Q £ eX p 0nen tial decay as t — > 00. By considering the Jordan canonical form 
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of A n , we obtain matrix polynomials Pj(t) such that e~ tAn = J2™ =1 Pj(t)e~ Kjt . Observe 
that F contains all the entries of e~ tAn B n B^e~ tA ™ and e~ tA ™C\aC n e~ tAn . The operator 
L x is an indefinite integral of e~ tAn B n B^e~ tA " while the operator Q% is an indefinite 
integral of e~ tA ™C\(jC n e~ tAn , hence L x and Q a x have entries in F ; moreover, the entries 
of (/ — L X Q^) _1 are quotients of determinants with elements in F . Hence by (9.32), 
logT n (2x) gives an element of F , so r n (x) itself is in a Liouvillian extension F of F . 

□ 

Theorem 10.2 Let ip be as in Theorem 9.2. 

(i) There exists a sequence of Unite rank matrices (A n )%L 1} with corresponding tau 
functions r n , such that ^logr n (2t) belongs to C(e~^ 1+1 / 2 ^, e~*) and r n (2t) ->■ r(2t) as 
n — > oo, uniformly on compact subsets of {t : $lt > 0}. 

(ii) Suppose further that k± is rational. Then there exists a positive integer N such 
that logr(2t) is periodic with period 2niN, and r n (2t) is given by elementary functions 
as in (10.4) below. 

Proof, (i) We introduce the finite rank matrices 

A n = diagonal [n t + 1/2, m + 3/2, . . . , k x + n + 1/2, 0, 0, . . .] (10.2) 

so that \\e~ tA - e ~ tAn || c i < e -(«i+«+i)H*/(l - e~ m ). Similarly, we cut down B and C and 
follow through the proof of Theorem 9.2 to produce the appropriate choice of W n (t,t) by 
the prescription of (9.8). Evidently, the eigenvalues of e~ tAn have the form e -~ t ( K i+ ( -+ 1 / 2 ) 
where I = 0, . . . , n. We observe that e -*("i+^+V2) belongs to C(e-^ 1+1 / 2 )*, e -*) for all £, 
so W n (t,t) likewise belongs to C(e-( Kl+1 / 2 )*, e - *). 

(ii) In this case, the set {m«i + m/2 + n£;m,n G Z;£ = 0,1,2,...} is a finitely 
generated subgroup of the rationals, and hence has a smallest positive element M/N, 
where M, N e N with M < N. Then all the terms N(ki + I + 1/2) are positive integers, 
so exp(-(t + 2nNi)A) = exp(-tA) for all 3fa > 0, hence r'(2t)/r(2t) is periodic. 

By Lemma 10.1, there exists a rational function r n such that 

^logr n (2t) = r n (e- t / JV ). (10.3) 
dt 

Suppose for simplicity that r n (z)/z has only simple poles; then from the partial fractions 
decomposition, there exist coefficients dj, fij and 7^ and bj, Cj such that 6 2 < Cj, real poles 
afc and a polynomial q n (z) such that (10.3) integrates to 

logr n (2t) = q n (e-*' N ) + £ a k log \ e -*' N - a k \ + £ & log(e~ 2t / N + 2b j e~ t / N + Cj ) 

k j 

+y: , i3 2 t^ e ~ t/N+ \ (10.4) 

3 y c 3 ~ b ] y c j ~ b j 
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When r n (z)/z has higher order poles, one likewise obtains expressions that are similar but 
more complicated. 



In view of Theorem 10.2(h) and the results of [8], it is natural to consider how the 
properties of r relate to integrability of Schrodinger's equation — /" + qf = Xf. 

Definition. Let q be meromorphic on C. As in [6,17], we say that q is algebro-geometric 
if there exists a non-zero R : C 2 — > C U {00} such that x i-> R(x; A) is meromorphic, 
A 1 — y R(x; A) is a polynomial, and 



Drach [6] observed that Schrodinger's equation is integrable by quadratures for all A 
only if q is algebro geometric. Conversely, if R(x; A) is as above, then 



gives a solution to Schrodinger's equation. In [6], Brezhnev catalogues several known 
special functions which give integrable forms of Schrodinger's equation. The theory extends 
to meromorphic potentials on compact Riemann surfaces by [13, p. 235] and [23, p. 1122]. 

The following result of Gesztesy and Weikard summarizes various sufficient conditions 
for a potential to be algebro-geometric. The initial hypothesis rules out variants of Bessel's 
equation. 

Theorem 10.3 [17] Suppose that —f" + qf = Xf has a meromorphic fundamental solution 
for each X and that either 

(i) q is rational and bounded at infinity; 

(ii) q is elliptic, that is, doubly periodic; or 

(in) q is periodic, with purely imaginary period, and q is bounded on {z G C : \$t,z\ > 
r} for some r > 0. 

Then q is algebro-geometric. 

We proceed to consider the cases (i),(h) and (hi) of Theorem 10.3, the linear systems 
(—A,B,C) that give rise to them, and the corresponding r functions. 

Proposition 10.4 Suppose that q satisfies Theorem 10.3(i). Then f has a rational Laplace 
transform and hence is the transfer function of a linear system (— A n , B, C) with a finite 
matrix A n . 



□ 



R'" + 4(q - X)R' + 2q'R = 0. 



(10.5) 




(10.6) 
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Proof. By a theorem of Halphen [17], the general solution of — /" + qf = Xf has the 
form fix) = Y^=i ( lj( x ) e ~ K: ' X i where qj(x) are polynomials. Hence there exist constants 
a/., not all zero, such that J2k=o a kf ( ' k H x ) = 0; so by taking the Laplace transform, and 
introducing the initial conditions, we can recover the rational function f(s); see [7, p. 15]. 
We recall that any proper rational function arises as the transfer function of a linear system 
that has a finite matrix A n , so /(A) = C(XI + A n )~ l B. 

□ 

Example (iii). In Theorem 10.2 (ii), q is a rational function of e~ f ^ N and under certain 
conditions gives rise to case (iii) of Theorem 10.3. In particular, q(t) = — 2sech 2 t is 
algebro-geometric, has period 2ni : is bounded on {z : > r} for all r > and arises 
from r(2t) = 1 + e~ 2t . This potential appears in the theory of solitons [4]. 

11. Realising linear systems for elliptic potentials 

Suppose that q is real, smooth and periodic with period one; introduce Hill's operator 
— -^2 + q{x) in L 2 (R). Then we introduce the Bloch spectrum, which is 

S B = {A G C : the general solution of - /" + qf = Xf has /GL°°(R;C)}. (11.1) 

One can show that when q is an algebro-geometric potential, Sb has only finitely many 
gaps; see [1,6, 24]. So we suppose that 

Sb = [A , Ai] U [A 2 , A 3 ] U . . . U [X 2g , oo) (11.2) 

with g gaps. The Xj are the points of the simple periodic spectrum, such that — /" + qf = 
Xjf has a unique solution, up to scalar multiples, that is one or two periodic. Let $ be 
the 2x2 fundamental solution matrix that satisfies 



1 

-X + q(x) 



$(x), $(0) 



1 
1 



(11.3) 



and let A(A) = trace $(1) be the discriminant of Hill's equation. We can characterize Sb 
as {A G R : A(A) 2 < 4}, and its components are known as the intervals of stability. 

The (hyper) elliptic curve C : y 2 = — YljLoiz ~ ^j) has genus g, and we can form the 
(hyper) elliptic function field S g = C(x)[y]. We therefore have a situation quite analogous 
to (2.5) and the Riemann surface £ of section 7. In this section, we consider the case of 
g = 1, where C is an elliptic curve which is parametrized by V . Hochstadt proved that 
g = 1 if and only if q(x) = c\ + 2V(x + c 2 ) where c\ and c 2 are constants; see [17]. 

Starting finite matrices, we formulate a version of the Gelfand-Levitan equation that 
is appropriate when <f>(x) = Ce~ xA B is periodic. (The equation (9.12) does not converge 
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when the functions are periodic.) A variant of this was used in [11] to solve the matrix 
nonlinear Schrodinger equation. Thus we will realise elliptic tau functions from linear 
systems. 

Definition. (Periodic linear system (— A, B, C; E)) Let A,B,C and E be finite square 
matrices of equal size; let e = ±1, and suppose that BC = e(AE + EA), BE = EB, 
EA = AE and exp27rA = /. Define 4>{x) = Ce~ xA B to be the scattering function for 
(—A, £?, C) and then introduce 

W(x, y) = Ce~ xA (I - e- xA Ee- xA Y X e- yA B. (11.4) 

We define the tau function to be 

r(x)=exp^y trace W(y, y) dyj (H-5) 

and let q(x) = —2^3 logr(x) be the potential function. 

We define r indirectly so as to accommodate the most significant applications. The 
definition retains the spirit of Theorem 9.2, on account of the following result. 

Lemma 11.1 (i) The matrices satisfy the Gelfand-Levitan equation 
-<p(x + y)+W(x,y)-e W(x, z)c/>(z + y) dz = W(x,y)E (0 < x < y < 2tt), (11.6) 

J X 



and 



-^-logdet(I - e- xA Ee~ xA ) = etrace W(x, x). (11.7) 
dx 



(ii) Let F be a differential field that contains all the entries of e~ xA . Then r(x) belongs 
to a Liouvillian extension of F, and r(x + 2n) = kt(x) where k = exp J Q 2?r trace W(y, y) dy. 
(in) Suppose moreover that e = 1 and 27r||</>|| 00 < 1. Then 

Proof, (i) One can check that 

/ e- zA BCe- zA dz = ee- xA Ee~ xA -eE (11.9) 

J X 

and it is then a simple matter to verify the integral equation (11.6). 
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By rearranging terms, one checks that 



trace W(x, x) = trace ((/ - e- xA Ee- xA )- 1 e- xA BCe~ xA " 

= e-^trace log(J - e~ xA Ee~ xA ) 
ax 

= e-^-logdet(J - e - xA Ee~ xA ). (11.10) 
ax 

(ii) By (i), r is given by exponential integrals of the entries of e~ xA . Note that W(x, y) 
is periodic in both x and y, so W(x, x) is periodic and hence J* traceVF(y, y) dy increases 
by the same amount as x increases through any interval of length 2tv. 

(iii) By repeatedly differentiating (11.6), and using periodicity, one derives the identity 

d 2 W d 2 W / d a „ „ dW 



+ 2 (^W(x, xj)<Kx + y) + W(x, O)0'(y) - -^-(z, O)0(y) 



fAr 2 <9y 2 V <7.r / ' (7// 

"27r ,q2tj; £)2wx fl2w fl2i 



r v (&w d 2 w\ , , a 2 w <9 2 w 



Since ABC - CSA = 0, we obtain 

O)0'(y) - -^-(x, O)0(y) = 0, (11.12) 

so (11.11) is a multiple of the original integral equation by — 2-^W(x, x). By the assump- 
tions on ||</>||oo, the solutions are unique, hence the differential equation is satisfied. 

□ 

By introducing infinite block matrices, we can extend the scope of Lemma 11.1. 
Clearly we can replace e in (11.6) by a diagonal matrix with blocks of ±1 entries on 
the diagonal. 

One can interpret the following result as saying that Lame's operator — -Ss + 2V 
has the scattering function proportional to sinx. Let ui± and ui2 be the periods, so that 
u = U2/0J1 has > 0; then let e± = V(ui/2), e<i = V((u>i + u>2)/2) and es = V(oj2/2); 
then let Jacobi's modulus be m 2 = (e2 — ez)f{ei — ez) and q = e luJ7T . To be specific, we 
choose wi = 2n and W2 = 2ni. Let A, B and C be the infinite block diagonal matrices 
with 2x2 diagonal blocks 

A = diagonal [ J] ~_ oo , C = A, (11.13) 

E = diagonalfq 2 !^]^, B = 2E. (11.14) 
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Proposition 11.2 (i) The functions 4>(x) = Ce xA B and W(x,y) of (11.4) satisfy the 
Gelfand-Levitan equation (11.6) and 



1 + q 2 

trace d>(x) = 4 - sin x; 

1 — q 2 



(11.15) 



(ii) The corresponding tau function is entire, belongs to a Liouvillian extension of the 
standard elliptic function field and satisfies 



2V(x) 



d 2 



dx 2 



logr(x). 



(11.16) 



Proof, (i) The matrices satisfy EB = BE, AE = EA and BC = AE + EA, so Lemma 
ll.l(i) applies. Note that the entries of E are summable, so E defines a trace class operator, 
hence the trace exists and a simple calculation gives (11.15). 
(ii) Observe also that 



det(/-q 2 l n| e- 2a;A ) = det 



1 — q 2 l n l cos 2x — q 2 l n l sin 2x 
q 2 Hsin2a: l-q 2 Hcos2x 



= l-2q 2|n| cos2x + q 



4|n| 



so one has 



det(J - e~ xA Ee~ xA ) = 4 sin 2 x J[ (1 - 2q 2n cos 2x + q 4n ) 2 ; 

n=l 

for comparison, by [25, p 135] the Jacobi elliptic function satisfies 

oo 

#1 (x) = 2q^ 4 sin x JJ (l - 2q 2n cos 2x + q 4n ) (1 - q 2 

n=l 

So we have an entire function 

t{x) = det(/ - e~ xA Ee~ xA ) = — 

Moreover, we have [25, p. 132] 



(11.17) 



(11.18) 



(11.19) 



9i{x) 2 

q 1/2 rL?=i(i-q 2n ) : 



(11.20) 



V(x) = -^\o g e 1 (x) + e 1 + ^log9 1 (x)\ x=1/2 , 



(11.21) 



hence we obtain (11.16). Let £ be the elliptic function field of functions of rational character 
on the complex torus C/(wiZ + w 2 Z). Then £ = C(V)[V], and by (11.21) £ has a 
Liouvillian extension Eq that contains 6\. 
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□ 

Theorem 11.3 Let r be an elliptic function. Then there exists a periodic linear system 
(—A, B, C; E), where A, B, C and E are infinite block diagonal matrices with 2x2 blocks, 
such that 

— log t(x) = trace W(x,x). (11.22) 

Proof. Any elliptic function is of rational character on C/(u\L + w 2 Z), and is the ratio 
of theta functions by [25, p 105], so 



Oix-dj) 

j=i 

where a\ + . . . + a m = b\ + . . . + b m . 

First we construct a periodic linear system with 9 as its tau function. For n = 0, let 
A = J/2, Eq = —iJ, B = il and Co = /, then (— Aq, B , Co; Eq) is a periodic linear 
system such that det(J — e~ xA ° E e~ xA °) = 2isina;. 

For n = 1, 2, ... , let A n = C n = J, E n = q 2n I and B n = 2E n ; then (-A n , B n , C n ; E n ) 
is a periodic linear system such that det(J — e~ xAn E n e~ xAn ) = 1 — 2q 2n cos2a; + q 4n . 
Hence we can introduce block diagonal matrices A = diagonal [Aq, A i, .. .] and E = 
diagonal [-Eo, E\, . . .], and so on to give a periodic linear system (—A, B, C; E) such that 

oo 

det( J - e- xA Ee~ xA ) = 2i sin x JJ (1 - 2q 2n cos 2x + q 4n ) 

n=l 

W{X) (11.24) 



"q 1/4 n^i(i-q 2n )' 

Next we replace (—A, B, C; E) by the terms (—A, e ajA B, Ce aj A ; e ajA Ee aj A ) which 
give Wj by (11.4); likewise we introduce (— A, e bj A B, — Ce bj A ; e bj A Ee bj A ) which give Wj 
by (11.4). We then form the block diagonal matrix 

®T=i {(-A) © {-A), e a > A B © e b > A B, Ce a > A © (-Ce b > A ); e a * A Ee a i A © e h * A Ee h * Ay ) (11.25) 
which gives the required W(x, y) = ®JL 1 Wj(x, y) © Wj(x, y) by (11.4), and we verify 

m 

trace W(x, x) = ^^trace Wj(x, x) + trace Wj (x, x)j 
j m 

= ^J2( l °S d ( x - a j)- l °sO(x-bj)) 

i=i 

= ^-\ogr{x). (11.26) 
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One can check that W satisfies (11.6) with e replaced by a diagonal matrix with diagonal 
entries ±1. 

□ 

12. Linear systems for potentials on hyperelliptic curves 

In this final section, we extend the analysis of section 11 to (11.2) in the case of g > 1. 
To obtain a model for the Riemann surface of C, we choose a two-sheeted cover of C with 
cuts along Sb, and introduce the canonical homology basis consisting of: 

• loops aj that start from [A2 3 ,oo), pass along the top sheet to [\2j-2, A2J-1], then 
return along the bottom sheet to the start on [A2 5 , 00); 

• loops fij that go around the intervals of stability [\2j-2, ^2j-i] that do not intersect 
with one another, for j = 1, . . . , g. 

Then as in [14, p 61], we form the g x 2g Riemann matrix [I; O] from the g x g matrix 
blocks 



x J 1 dx ] 9 

y J j,fc=i' 



and O 



1 dx 



y J j,fc=i 



Then the corresponding Riemann theta function is 

0(s I O) = exp(z7r(On, n) + 2ivi(s, n)). 



(12.1) 



(12.2) 



Example. Suppose that g = 2 and let 



n 



a b 
b d 



(12.3) 



where > 0, Qd > and b G Q. Then choose p G N such that pb G Z. One can easily 
check that 

p-i 

0(3, t\Q)= e^+^+^V^V^fcs + r I p 2 a)%t + I P 2 d). (12.4) 

r,/j=0 

Proposition 12.1 Suppose that q is a periodic potential with g spectral gaps, as above, 
and that 0( | Q) is a finite sum of products of Jacobi elliptic functions. Then there exist 
N < 00, Xj G R, aj G C with Scr,,- > 0; and block diagonal matrices (A,-, Bj, Cj; Ej) with 
2x2 diagonal blocks for j = 1, . . . , N, such that 9(x — Xj \ <jj) is the tau function of 
{—Aj, Bj, Cj; Ej) and q belongs to the field C(6(x — Xj \ <Jj);j = 1, . . . , N). 

Proof. We introduce the product of real ovals 



{ g ( A (^) + V /4 - A (^) 2 ) 5 =1 : A 2j-i < ^ < A 2j : J = 1, • • • , 9} 



(12.5) 
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which has dimension g. McKean and van Moerbeke [24, p260] considered the manifold M. 
of all the smooth real one-periodic potentials such that the corresponding Hill's operator 
has simple spectrum {Ai, . . . , A2 5 }. Using the KdV hierarchy, they introduced a 1 to 
1 differentiate map from M onto T s , where the tangent vectors on M are differential 
operators. Let A be the lattice generated by the columns of [I; O], and note that C 9 / A is 
the Jacobi variety of C. They showed that q extends to an abelian function on C 9 which 
is periodic with respect to A, hence gives a function of rational character on C 9 / A. The 
extended function q belongs to £ g , hence is a theta quotient. Moreover, translation on the 
potential is equivalent to linear motion on C 9 at constant velocity. 

Thus they solved the inverse spectral problem explicitly, by showing on [24, p. 262] 

that 

A 0(X - u*/2 | fl)Q(X - ou**/2 | ft) 

= h £j Q ( x - ^ 2 1 n W x - w ~/2 1 n) (12 ' 6) 

where X = (x\, . . . , x g -i, ax + b) has fixed, while a; varies, and the con- 

stants Sj, LOj, a;**, and are notionally computable. 

By hypothesis, each factor Q(X — u* /2 \ ft) may be written a a finite sum of products 
of functions such as 9(ax + Cj \ dj), and we can apply Theorem 11.3 to each such factor. 

□ 

Weierstrass and Poincare developed a systematic reduction procedure for such elliptic 
functions of higher genus, so we can describe the scope of Proposition 12.1. The Siegel 
upper half-space is 

S g = {ft G M gXg (C) : ft = Q*;QQ > 0}. (12.7) 



Let X and J be the 2g x 2g rational block matrices 



X 



a (3 
7 5 



J 



-I 

1 



(12.8) 



such that XJX 1 = J; the set of all such X is the symplectic group Sp(2g; Q). Now X is 
associated with the transformation <fx of S g given by 

cpx{ft) = (al + l3fl)- 1 ( 1 l + 8fl), (12.9) 

thus Sp(2g; Q) acts on S g . 

Proposition 12.2 [1] (i) Suppose that fl can be reduced to a diagonal matrix by the 
action of the symplectic group. Then 0( | ft) can be expressed as a sum of products of 
Jacobian elliptic theta functions. 
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(ii) Condition (i) is equivalent to C being a N -sheeted covering of the one-dimensional 
complex torus for some N. 

(in) The orbit of Sp(2g, Q) that contains ii is dense in S g . 
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